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Abstract 


In this thesis the evaluation of convolution integral using continuous wavelet transform 
has been studied for two time domain signals The results of the convolution product 
thus obtained has been compared with that of the direct convolution integral method 
The results indicate that the choice of the wavelet basis function and the corresponding 
dilation and translation parameters play an important role m the correct evaluation of 
the convolution integral 
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Chapter 1 


Introduction 


In recent times wavelet analysis (also called wavelet theory or just wavelet) has at 
tracted much attention in signal processing It has been successfully applied in many 
applications such as transient signal analysis, image analysis communication systems, 
and other signal processing applications Its attraction is due to its ability to anal 
yse rapidly changing transient signals It can be shown that every application using 
the Fourier Transform (or FFT) can be formulated using wavelets [1] to provide more 
localized temporal (or spatial) and frequency information 

In signal processing wavelets are very useful for processing non stationary signals 
When analysing such signals with Fourier Transform any abrupt change m time in the 
signal x(t) IS spread out over the entire frequency axis in X(f) where X(f) is the Fourier 
Transform of x(t) Therefore analysis adapted to the non stationary signals requires 
methods other than the Fouiier Transform in order to know the frequency content 
of the signal with change in time that is we need a two dimensional time frequency 
representation of the signal x(t) Short Time Frequency Transform(STFT) is one such 
lepresentation which uses time window function of constant width But this has some 
drawbacks as far as the frequency and time resolution is considered [1] 

In the STFT analysis of a non stationary signal there are two conflicting require 
ments the window width must be long enough to give the desired frequency resolution 
but must also be short enough so as not to blur the time dependent events Given a 
window function g(t) and its Fourier transform G(f) the frequency resolution A/ and 
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2 


the time resolution Aj are defined as [2] 

’ S\au)?df 

and 


(1 1 ) 


( 12 ) 


The frequency resolution and the time resolution are related by the uncertainty 
principle as 


AtA/>^ (13) 

471 

1 e their product is lower bounded Equationl 3 indicates that one can only trade 
frequency resolution for time resolution or vice versa More important is that once a 
window has been chosen for the STFT then the time-frequency resolution given by 
equation 1 1 and equation 1 2 is fixed over the entire time frequency plane 

The Wavelet Transform representation for non stationary signals overcomes this 
limitation through the use of varying window widths by choosing different scales for 
baby wavelets [3] Let ipij:) be the mother wavelet which has Aj = Sf and At = 5* then 
ipit/a) baby wavelet will have Af = 5f fa and At — a5t where a is the scale parameter 
By varying a various time and frequency resolutions can be obtained although of 
course not both simultaneously 

The wavelet theory has been applied in various fields In computer vision multires- 
olution techniques have been used for various problems ranging from motion estimation 
to object recognition [4] Images aie successfully approximated starting from coaise 
vcision and going to a fine lesolution \eisioii 

One of the most impoitant ipplications for wavelets is in signal compression /cod- 
ing of image and speech Applications of wavelet decomposition in numerical analy- 
sis e g for solving partial differential equations seem very promising because of the 
zooming property which allows a very good representation of discontinuities unlike 
the Fourier Transform [4] Statistical signal processing using wavelets is another appli- 
cation of wavelets 
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In signal compression the multiresolution technique resolves the input signal f{t) 
into components of different resolutions A. coding scheme can then asign more bits 
to represent the lower detail sequences which have shorter lengths due to successive 
subsamplings and less to the higher detail sequences in order to achieve compression 

One particular application of continuous wavelet transform is to perform linear hi 
tering As we will see in chapter 3 it is possible to obtain convolution of two signals 
with the continuous wavelet transform The complexity of this algorithm is proper 
tional to N If however we achieve the convolution with the FFT algorithm the 
complexity of the algorithm is Nlog 2 N the computing time is thus longer than the 
one obtained with the Wavelet Transform if we concentrate the energy in very few 
coefficients 

In the present thesis work the computation of convolution using continuous wavelet 
transform has been undertaken and the results have been discussed 

1 1 Organisation of the thesis 

The thesis has been presented m four chapters In chapter 2 the wavelet theory in 
general and different wavelet transforms alongwith multiresolution analysis has been 
discussed Chapter 3 presents the theory and the algorithm for computing convolu 
tiuon using continuous wavelet transform The results and the conclusions are given 
in chapter 4 



Chapter 2 


Wavelet Theory 


The mathematics associated with building a model for a signal system or process 
with a set of special signals ’ or ” wavelets as known popularly is called as the 
wavelet theory The wavelets must be oscillatory and have amplitudes which quickly 
decay to zero in both positive and negative directions Wavelet theory represents 
things by breaking them down into many interrelated pieces when the pieces are 
scaled and translated wavelets this breaking down process is termed as the wavelet 
decomposition or wavelet transform Wavelet reconstruction or inverse wavelet 
transform involves putting the wavelet pieces back together to retrieve the original 
object or process Wavelet theory involves the scaling or warping operation and the 
translation or shifting operation If the translation operation is combined with the 
scaling action then the composite operation is referred to as an affine operation [5] 
The ifhne operation simultaneously scales and translates the independent variable 
M ithematically the affine operation maps the original independent variable t into a 
new independent variable t by the formula 


t' = 


t — b 


a 

where a is the scale and b is the translation 


(2 1) 
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2 1 Wavelet Transform 

Like Fourier analysis, wavelet analysis deals with the expansion of functions in terms of 
a set of basis functions but unlike Fourier analysis wavelet analysis expands functions 
not in terms of the trigonometric polynomials but in the form of translations and 
dilations of mother wavelet The Fourier transform of a given signal gives the frequency 
content of the signal However it fails to provide information about the time evolution 
of the frequency content i e , how the frequency contents of the given signal change with 
time This time frequency localization is achieved by the Wavelet Trans form{WT) 
Thus for the analysis of non— stationaary stganl the wavelet transform provides an 
alternative to the classical Short Time Fourier Transform (STFT) or Gabor transform 
In contrast to the STFT which uses a single analysis window the WT uses short 
windows at high frequencies and long windows at low frequencies, so that the time 
resolution becomes arbitrarily good at high frequencies while the frequency resolution 
becomes arbitrarily good at low frequencies For example two very close short bursts 
can always be eventually separated in the analysis by going up to higher analysis 
frequencies in order to increase time resolution This kind of analysis of course woiks 
best if the signal is composed of high frequency components of short duration plus low 
frequency components of long duration which is often the case with signals encountered 
in practice There are four different types of wavelet transform [3] 

2 11 Continuous Wavelet Transform 

The Continuous Wavelet Transform [CWT) of a signal f{t) maps the one 
dimensional function into a two dimensional function of scale a and translation b and 
IS defined as, 

CWT{a b) = ^j f{t)^ (^) 
where the signal f{t) is assumed to be square integrable denoted as 
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fit) e V{R) (2 3) 

meaning 

j P{t)dt < oo (2 4) 

and 0(t) IS the basic (or mother) wavelet (t) denotes the complex conjugate 
of ipit) ip{t) can be a real or a complex valued function which satisfies the following 
conditions 


and 



1 7/'(i) dt < oo 


(2 5) 




I ^(w) P 
I w I 


dw < oo 


(2 6 ) 


Where ^^(’w) is the Fourier transform of ipit) The first condition implies finite 
energy of the function t/i(i) and the second condition is the admissibility condition 
In the CWT the variable t scale a and shift b are all continuous Fig 2 1 shows the 
schematic of obtaining the wavelet domain coefficients which are computed for each 
particular scale and translation (a b) value Thus we see that the continuous wavelet 
transform is an analysis filter which breaks a signal down into component pieces and 
these pieces art examined or operated instead of the original function 

The mver se continuous wavelet transform maps a surface in the scale translation 
plane H(a b) into the one dimensional time domain to reconstruct the original signal 
/(f) ind IS defined as 



da db 
a2 


(2 7) 


Wavelet Transform 
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The inverse transform thus creates the original function by summing appropriately 
weighted scaled and translated versions of the mother wavelet The weights are 
the wavelet coefficients CWT{a 6) Figure 2 2 illustrates the reconstruction of original 
signal f[t) from the wavelet domain coefficients The inverse transform or the recon 
struction filter is a synthesis filter The synthesis filter puts the pieces back together 
again A characteristic of the inverse transform is that it is non unique in the sense 
that several different wavelet transform domain representations can be inverted with 
respect to the same mother wavelet to create the time function f{t) That is the 
inverse transform may be a many to one operator This non-uniqueness is due to the 
increase in the dimensionality 

2 12 Discrete Parameter Wavelet Transform 

In the continuous wavelet transform both scale a and translation b are continuous van 
ables and there is a redundancy in the CWT representation of (a b) The Discrete 
Parameter Wavelet Trans form{DPWT) of a continuous signal takes into consider 
ition this fact and it is defined as [3] 


DPWT{mn)—aQ J f{t)i}{aQ'^t — nbo)dt (2 8 ) 

where the parameters a, and b are discretized to a = and b = nboa'Q with uq bo 
the sampling intervals and m n mtegeis Both f{t) and 'tp(aQ^t) are still continuous 
For computational efficiency ao - 2 and bo = I are commonly used resulting in a 
bin uy dilation of 2~ "-and a dyadic tianslation of 2"‘n Equ ition 2 8 can also be written 
in the inner product form as 

DPWT(m n) =< f{t) ibr^nit) > (2 9) 

m 

where '0mn(t) = “o ^ 


Inverse Wavelet Transform 
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Figure 2 2 Inverse Wavelet Transform Structure 
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The original time signal can be reconstructed from the DPWT{m n) as 

fit) = c E E DPWT{m (2 10) 

m n 

where c is some constant dependent only on '0(t) In case of oversampling the above 
equation will not only hold but moreover it is possible to have a non unique represen 
tation of f(t) with respect to the same That is different sets of (ao bo) resulting 
in different iJranii) will produce different DPWT{m n), each capable of recovering the 
signal /(t) In contrast if the sampling is sparse it may not be possible to recover /(t) 
at all 

2 13 Discrete Time Wavelet Transform 

The Discrete Time Wavelet Transform{DTWT) is a time discretization of DPWT 
with t = kT and the sampling interval T = 1 It is defined as 

DTWT{m n) = Uq ^ E (2 H) 

k 

We note that for ao = 2 there is an output only at every 2”^ samples when 2~'^k is 
an integer 

2 14 Discrete Wavelet Transform 

The discrete wavelet transform is defined as 


DWT{m n) = 2-t E s{k)iP{2-'^k - n) (2 12) 

k 

where the discrete wavelet i/}{k) can be but not necessarily a sampled veision 
of 1 continuous counterpart When i){k) is a discretization of a ip{t) the D\V T is 
identical to the DTWT 
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2 2 Multire°olution Analysis 

Multiresolution analysis is the decomposition of a signal f{t) into components of dif 
ferent scales (frequency) of 2~’^ m E Z [6] Associated with each scale is a subspace 
Vm These subspaces are time functions which satisfy the following conditions stated 
below [4] 

i Containment 

{0} ^2 c Ml C Mo C M _1 C M -2 -4 L2{R) (2 13) 

ii Upward Completeness 

Uv;n=L2(R) (2 14) 

m&Z 

III Downward Completeness 

n = { 0 } (2 15) 

IV Scale Invariance 


fit) G Mn « fi2^t) e Mo 


(2 16) 


V Shift Invariance 


f{t) eVo<=> fit - n) G Mo for all n e Z (2 17) 

VI Existence of orthonormal scaling functions 

There exists a scaling function ^{t) e Vq such that the set 

{<A„„(t) = 2'¥0(2-"'i-r.) 716 2} (2 18) 

is an orthonormal basis that spans Mn 

Figure 2 3 shows the MRA decomposition scheme and figure 2 4 gives the structure 
for reconstruction The multiresolution analysis scheme decomposes fit) into the de 
tail components m = 1,2 L and the coarsest approximation via stages of 
identical low pass filters gil) and high pass filters /i(/) and subsamphng^ the outputs 


mm 

Ml 




I } 


iiUH 
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Figure 2 4 MRA Reconstruction Scheme 
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The decomposition continues until it reaches conceivably the dc component 
of f{t) The subsamplmg by two produces a time compression by two and identical 
'g{l){h{l)) will have bandwidths of the lower(upper) half of the previous g{l){h{l)) The 
details aie the DPWT coefficients of f(t) le — DPWT{m n) under certain 
conditions[3] The choices of g{l) and h{l) determine (f)(t) and The filters g[l), 
h{l) g{l) h{l) are related to scaling function and wavelet function as [3] 



(2 19) 

1=0 


P-1 


ip{t) ^ ~ 1) 

1=0 

(2 20) 

h{l) = {-iyg{p - 1 - 1) 

(2 21) 

h{l) = h{p — 1 — 1) 

(2 22) 

W) = - 1 - 0 

(2 23) 


where the delay p - 1 is the filter order The above- scheme is the M allot algorithm 
foi MRA also known as the Fast Wavelet Transform when attached with a WT 
interpretation 

The algorithm has complexity 0{N) and is therefore faster than FFT The com- 
plexity IS independent of the number of stages of decomposition Since a subsamplmg 
by two follows each stage the total no of computation is 

T = To + y + ^+ =2To (224) 

whcie To IS the no of computations for the first stage 


Chapter 3 


Convolution And Wavelet Transform 


Amongst several applications of wavelet theory we consider the problem of computing 
the convolution product of two time signals using wavelet transform 

3 1 Convolution Integral 

The convolution integral of two functions f{x) and g{x) is defined as 

/ OO 

f{u)g{x — u)du (3 1) 

■OO 

and is denoted as 

h{x) = f{x) *g{x) (3 2) 

For Linear Time Inuaiiant {LTI) systems the convolution integral represents the 
output h{x) of the system in terms of the input signal f{x) iiid the system impulse 
lesponse g{E) [7] Convolution integral is the basic time domain formula in evaluating 
the system response 

3 2 Convolution Theorem 

If Fi{w) is the Fourier Transform of fi{x) and F^iw) is the Fourier Transform of f 2 {x) 
1 e 


/i(s) — 5- Fi{w) 


(3 3) 
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and 


f2(x) ^ F2 {w) (3 4 ) 

then 

y(x} = fi(x) * f 2 {x) < — )■ Fi{w)F2{w) (3 5 ) 

1 e 

Y{w) Fi(w)F2{w) (3 6 ) 

where y{x) is the convolution product of fi{x) and f2{x) and Y{w) is the Fourier Trans 
form oiy{x) Convolution Theorem is an important theorem applied in evaluating the 
system response in the frequency domain 

3 3 Convolution Using Continuous Wavelet Transform 

We will examine here the computation of a convolution using continuous wavelet trans 
form as presented in [8] Let us consider the convolution product of two functions f{x) 
ind g{x) 


/ OO 

f{u)g{x — u)du 

-OO 

Iiitioducing two real wavelet functions 0(x) and \;(/) such that 


(3 7 ) 


C=r (3 8 ) 

Jo u 

IS defined Here 0(r/) and x{u) denotes the Fourier Transforms of the wavelets '4>{x) 
and x{x) respectively 

The Wavelet Transform of g{x) with respect to the mother wavelet 'i/)(x) is 


W,{a b) - 



(3 9 ) 
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The inverse wavelet transform of Wg{a b) with respect to the wavelet function xi^) 
can be written as 


1 roQ fOO 1 ( r — 

Substituting g(x) in equation 3 7 we get 


X — b\ da db 


/ oo 1 rco roo 1 

fi^)r L / -^Wg{a,b)x 

-oo O Jo J —CO \j a 


' X — u — b\ da db 


du 


a I a^ 


(3 10) 


(3 11) 


h{x) = ^f Wg{a b)db f fiu)xl- — - — cfw (3 12 ) 

L' Jo Cl J—oo J—oo \ a / 


X — u — b^ 


We define 


Xix) = xi-x) 

so that the WT of f{x) with respect to the wavelet function x{x) is 


Now the integral from equation 3 12 can be written as 


/ OO 

f{y)x 

-oo 


^ L ' 

X — U -- 0 


du 


/ OO 

fiu): 

-oo 


U 


-{x-b) 


du 


/•OO 

= / f{u)x 

u — 

{x-b)' 

du 


J—oo 

, 

Cl 



ix 

II 

1 

1 

-5)1 

du 

J—oo 


L “ 

_ 

. 


(3 13) 


(3 14) 


\/5 


I e 


/ C50 

f{u)x 

-CO 


' T — u — b' 


du = ^/aWf(a, x-b) 


(3 15) 


Using the above expression we may write from equation 3 12 

1 r°° da /■'» 


h{x) = ^ 

C Jo a^ J-oo 


(316) 
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but the second integral of equation 3 16 can be represented as the convolution product 
of Wf{a x) and Wg{a x) Then we get the final result 

1 ~ fin 

h{x) = - / Wf{a x) ® Wg{a x)^ (3 17) 

L/ JO 

Therefore in order to compute a convolution with the continuous wavelet transform 

• we compute the wavelet transform Wf(a b) of the function f{x) with the 
wavelet function x(x) 

• we compute the wavelet transform Wg(a b) of the function g{x) with the wavelet 
function ■ 0 ( 2 ;) 

• we compute the convolution product of the wavelet transforms Wf(a x) and 
Wg{a,x) 

• we sum this convolution for all scale values from 0 to co 

The Wavelet Transform permits us to perform any linear filtering Its efficiency 
depends on the number of terms in the WT associated with g{x) for a given signal f{x) 
If we have a filter where the number of significant coefficients is small for each scale 
the complexity of the algorithm is proportional to N For a classical convolution the 
complexity is also proportional to N but the number of operations is also proportional 
to the length of the convolution mask The mam advantage of the present technique 
lies in the possibility of having a filter with long scale terms without computing the 
convolution on r large window If we achieve the coniolution with the FFT algorithm 
the complexity is in Nlog 2 N the computing time is longer than the one obtained with 
the wavelet transfoim if we concentrate the energy on very few coefficients 

3 4 Convolution Using DPWT 

In case of Discrete Parameter Wavelet Transform {DPWT) the wavelet transform 
IS not taken for all scale and shift values but computed only for fixed number of these 
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values ds decided by the sampling grid The convolution of two continuous signals 
using discrete parameter wavelet transform has been derived and presented below 
Considering again two real wavelets and x{^)> the DPWT of g(x) with respect 
to the wavelet function ijj^x) is given by 

/ OO 

g{x)'lpmn{x) (3 13 ) 

-OO 

where iJmnix) = 2-T'tlj{2-^x - n) for dyadic sampling We can recover g(x) with the 
wavelet function xi^) so that 

g{x) ^ c'£^Y^Wg{m n)xmnix) (3 19 ) 

m n 

where Xmn{‘^) — 2 — x(2"’^a:-n) for dyadic sampling and c is some constant dependent 
only on the wavelet x(^) 

Substituting g(x) from equation 3 19 in equation 3 7 we get, 

/ OO 

f{y)Xmn{^ - u)du (3 20) 

"OO 

jri ii> 

Denoting x(x) = x(-3:), the DPWT of f{x) with respect to the wavelet function 

Xix) IS 


Wf{mn)= / fix)Xmn{'^)dx (3 21) 

J —OO 

where Xmn{-^) = 2“^y(2~"'c - n) Now consider the integral from equation 3 20 

fiu)Xmn{x - u)du = J^oo f{n)2-TXmn - u) - u] du 

= 2-t f{u)Xmn [- {n - 2 - " - y-)}] dy 


- 2-t f[u)x [n - 2-"^ (x - u)] du 


= 2 -^ !Zo fiy)x - n)] du 
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1 e 

/ OO 

/ (u)Xmn(x - u)du = Wf (m 2~'^x - nj (3 22) 

From equation 3 20 and equation 3 22 we get the final expression for h(x) as 

h{x) — n)Wj{m 2~^x - n) (3 23) 

m n 

The above equation is similar to that obtained with the continuous wavelet transform 
And we can obtain the convolution product of two signals using the DPWT 


Chapter 4 


Implementation And Results 

4 1 Choice of Wavelet Basis 

In implementing the computation of comolution through i\a\elet transform we need 
two leal wavelets (/ (t) and x(0 such that C — defined In the present 

method we have chosen and Y(t) such that 

ip{t) = x{t) (41) 

where ip{t) is an e\cn function of time The wavelet ?/)(f) chosen m our implementation 
is the second den\ative of Gaussian wavelet which is an even function of time and it 
is defined as [3] 

= (1 - (4 2) 

md its F T is gupii by 

4'(u;) = x/Trui^e (13) 

These functions ire plotted in the hguie 4 1 uid hgiire 4 2 
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4 2 Implementation of Convolution Integral 

From equation 3 9, the continuous wavelet transform of the signal g(t) with respect to 
the wavelet function ip{t) is given as 

W,ia b) -L r (—] dt (4 4) 

yj d J — oo y CL J 

substituting from equation 4 2m the above expression we may write 

Wg{a b) ^ f g{t) 1 - (- — (4 5) 

CL 7—00 y u / 

and similarly the continuous wavelet transform of the signal f{t) with respect to 
the wavelet function Lp{t) is given as, 

Wf{ab) = -j= /(i) 1- e-^-r-dt (4 6) 

sy CL 7—00 y <1 y 

Now changing the time variable t = uin the equation 4 6 we get 

Wfia,b) = ^f f{u) 1-f^^ — e~^"~^du (4 7) 

y CL J —oo \ ^ J 

substituting equation 4 5 and equation 4 7 in equation 3 16 we get the expression 
for h{x) as 
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Equation 4 9 involving four integrals was evaluated using Gaussian Legendre m 
point formula method [9] Two different sets of input signals f{t) & g(t) were taken 
The first set was simulated with a ramp signal and a gate signal and the other set with 
two exponentials having different time constants 

4 2 1 Ramp and Gate Signals 

The functions f{t) and g{u) were taken to be gate signals and ramp signals respectively 
as defined below 


fit) - 1 0 < i < 5 
= 0 elsewhere 


md 


(410) 


giu) = u 0 < tt < 10 
= 0 elsewhere 


(411) 


Therefore in equation 4 9 substituting equation 4 10 and equation 4 11 and chang- 
ing the limits accordingly we get 


j•oo TOO rS rW \ 

/l( l) = 7^ / / / / "7'* 

C Jo J—ooJo Jo O, 


t b — X 


(^) 


'u — h'' 




du dt db da 


(4 12) 


We have calculated the convolution integral of a gate signal as given by equation 
4 10 and a ramp signal as given by equation 4 11 by using direct integration, using 
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time domain convolution method graphically Considering the appropriate limits for 
diffeient overlapping regions of the functions f(u) and g(x u), the convolution product 
of f(x) and g(x) can be obtained as [7] 

h{x) = 0 a: < 0 

= 0 < X < 5 

= 5x - 12 5 , 5 < X < 10 (4 13) 

= + 5x + 37 5 , 10 < X < 15 

= 0 15 <x 

The convolution integral so calculated is shown in the figure 4 3 

4 2 2 Exponential Signals 

Two exponentials having different decay constants were also taken for performing con 
volution using continuous wavelet transform The functions f{t) and g{t) are defined as 

f{t) = exp{-at)u{t) (4 14) 

md 


g{t) - exp{-pt)u{t) (4 15) 

Therefore in equation 4 9 substituting equation 4 14 and equation 4 15 we get 


■) fOD / CO fOO roo 1 

h{i) = ~ / / / ^exp{-at)exp{-pu) 

C Jo J '-oo J —oo J —oo CL^ 


't + b 


^ t+t-j ^ 



(4 16) 


du dt db da 



Figure 4 4 Convolution of Ramp and Gate Signal Using Wavelet Transform 
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Figure 4 6 Convolution of Ramp and Gate Signal Using Wavelet Transform 











key {gauss_amin_amax_bmin_bmax) 
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Figuie 4 8 Convolution of Ramp and Gate Signal Using Wavelet Transform 
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Figure 4 11 Convolution of Ramp and Gate Signal Using Wavelet Transform 
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The time domain, convolution of the exponential signals given by equatin 4 14 and 
equation 4 15 has been calculated and the expression for the convolved output is 
given by 

h(.a:) {^~a) ^ ^ P 

= X exp{-ax)u{x) if a = 13 (4 17) 

The direct convolution of these signals obtained eith the above equation is shown 
m figure 4 12, and figure 4 13 

4 3 Results And Discussion 

We have calculated the direct convolution integral of a ramp signal and a gate signal 
graphically and is given by the equation 4 13 Figure 4 3 shows this as a function of 
the variable x The direct convolution of the two exponentials, with different decay 
constants a and P as given m equation 4 14 and equation 4 15 has also been calculated, 
and the closed form expression for the same is given in equation 4 17 The figure 4 12 
and figure 4 13 shows the convolved output for different value of P 
The convolution integral so calculated for both sets of input signals has been compared 
with that obtained using continuous wavelet transform (equation 4 12 and equation 
4 16) The comparison is given below 

1 The parameters to be chosen for the calculation of the integral given by equation 
4 12 aie the lowei and upper limits of scale o and tianslation b Those for t and 
a ire known since we have used known signals 

2 An inspection of equation 4 12 shows that the integrand is always finite including 
a = 0 rnd a = oo for well behaved signals f{t) and g{u) The rise in the integrand 
rt a = 0 because of the of term in the denominator is compensated by the 
exponential term in the numerator, in which a appears in the denominator of the 
argument (which is negative) of the exponential function Thus in the integration 
the lower limit of a can be set as zero without any problem Large values of a are 
unlikely to contribute much because of the term in the denominator Small 
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values of a also are unlikely to contribute much because of the structure of the 
exponential term m the numerator However fixing of the upper cutoff (nmax) 
needs some experimentation These results are discussed later in this section 

3 The upper and lower limits of the scale b vary with the value of the scale a We 
have generated the integrand of the inner integral (integrand over the variable h) 
m equation 3 16, at different discrete values of a and x and find that the range 
of b over which the integrand assumes significant values vary with both a and x 
A range of 50 to +50 would be adequate over the range studied (upto a = 10, 
for different values of x) 

4 A severe restriction in our calculation has been the number of points that we use 
for evaluating the final integral (equation 412) The program takes a lot of time 
if the number of points are increased beyond 20 or at best 30 Thus increasing the 
range of b to 50 to +50 is not necessarily more accurate than a smaller range of 
b, because as we increase the range we consider a more complete range of values 
of the integrand, but since we do not take proportionally large number of points, 
discretization error enters 

5 Figure 4 4 to figure 4 6 shows some of the better results we have obtained Figure 
4 4 shows a graph with limits for scale a from 0 to 100 and for the translation 
b from —20 to +20 It shows a maximum value of approximately 35 at a: « 10 
very close to that in figure 4 3 There is a little bit of ripple and the graph starts 
a little below zero at a; = 0 but the overall agreement with figure 4 3 is very 
good Figure 4 5 is also in good agreement with figure 4 3 , but the agreement is 
not as good as it is between figure 4 3 and figure 4 4 The value of „ m figure 
4 5 IS 0 4 and not zero as in figure 4 4 The small yet noticeable deteiioration 
in agreement shows that the contribution of n = 0 to 0 4 is not to be ignored 
Figure 4 6 shows a graph that has same parameters as figure 4 4 but the range 
of b IS from 50 to +50 The agreement with figure 4 3 is not as good as it is 
with figure 4 4 A wing appears on the right and the maximum amplitude is 


33 


(lepiessed The lesult shows that if number of points for integration does not 
increase propoitionally discretization error increases as lange increases The 
deterioration however is not drastic Increasing the range of 6 to 100 to +100 
and cleci easing the integration points to 10 from 20 deteriorates the agreement 
We note that inclusion of the lower limit of the scale a from = 0 0 to 0 4 
improves the lesult The tact that even with amjn = 0 4 we get very satisfactory 
lesult agiees with the conclusion drawn earlier that lower values of a do not 
contribute significantly We had anticipated that an upper limit of the scale a 
= 10 would be adequate because of the factor in the denominator of 
equation 4 12 But our results indicate that Umax = 100 is necessary to obtain 
good ht 

6 If we insist on a smaller amax(w*-2 10) in order to avoid discretization error we 
obtain t highly oscillatoiy function h{x) with a maximum at around a = 10 
but the oscillations mask the actual value of h{x) Such a graph is shown in Fig 
4 7 to Fig 4 8 The maximum shoots off to a very large values, much larger than 
the actual value of 37 5 which is nearly achieved m Fig 4 4 when Omax is set at 
100 

7 We note th it the oscillatory nature of calculated h{x) obtained when amm = 0 
md n„,ax = 10 (figure 4 7 and 4 8) disappeais when Um n = 0 3 and above (figure 
1 9) This combined with the fact that oscillations disappears with Cmm = 0 if 
<''ni IX = 100 sugg( Sts th it the values of a between 0 0 md 0 3 or 0 4 aie responsible 
foi the oscillitoiy behivioui of h{i) In both cases with am n = 0 and a„ia = 100 
md (imin = D i 111<-1 “mLx = 10 sincc the nuinbei of points for evaluating the 
mtegi il IS birae sin ill vilues of a ue ivoided 

S We further note that m Fig 18 if wc draw a mean curve by avenging the 
oscillitions w( g(t ireisonable ipproximation to Fig 4 3 The maximum value 
o((uis IS It should be it around j; = 10 but its magnitude is somewhat less 
th in that in Fig t 3 ( '^ 30 instead of 37 5 ) This is not true if Omm ~ 0 
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9 We have also computed h{x) for amax = 20 with amm = 0 and other small values 
VIZ amm = 0 1 to 0 5 and find highly oscllatory nature m h{x) when Umm = 01 
or 0 2 but the oscillations disappear at = 0 4 or 0 5 

10 The conclusion that smaller values of a are responsible for the oscillatory nature 
of h{x) is supported by the result that increasing the number of points used in 
the integration for fixed values of other parameters increase oscillations Fig 
4 10 and 4 11 demonstrates this effect 

11 Convolution integral of two decaying exponentials the decay constant of the first 
being 1 and that of the second being 1 10 20 and 50 are calculated by method of 
continuous wavelet transform with the limits for scale a from 0 to 50 and that 
for the translation b from 20 to +20 The direct convolution using Eq 4 17 
has been shown in the Fig 4 12 and 4 13 The convolution using continuous 
wavelet transform (Eq 4 16) has been shown in the Fig 4 14 to Fig 4 17 
The agreement of these plots with Fig 4 12 and 4 13 is quite good Detailed 
investigations on choice of the limits for the scale a and the translation 6 are 
however not done 

4 4 Conclusion And Future Work 

Convolution integial of two signals can be satisfactorily calculated using the method 
of continuous wavelet transform with suitable choice of the wavelet function and with 
appropriate choice of the range of dilation and translation parameters This has been 
illustrated by evaluating the convolution of a ramp and a gate signal and of two ex 
ponentials, with the second derivative of Gaussian wavelet Spurious oscillations in 
calculating the function h{x) are observed under certain conditions which suggest 
that small values of the scale parameter a ( between 0 to 0 2 ) cause the oscillatory 
behaviour The reason for these oscillations are, however, not fully clear Further in 
vestigations of this problem for different choice of the range of parameters a and b can 
be taken up as the future work 
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Figure 4 16 Convolution of Two Exponentials with Different q;(/? = 1) Using Wavelet 
Transform 



Figure 4 17 Convolution of Two Exponentials with Differnet q;(/ 3 = 1) Using Wavelet 
Transform 
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